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EFFECT OF COULOMB SPLINE ON ROTOR DYNAMIC RESPONSE 
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A rigid rotor system coupled by a coulomb spline is modelled and 
analyzed by approximate analytical and numerical-analytical methods. 
Expressions are derived for the variables of the resulting limit cycle 
and are shown to be quite accurate for a small departure from isotropy. 

SYMBOLS 

a,b semi-major, minor axis of the elliptical response 
a m (a+b)/2 

a d (a-b)/2 

c external viscous damping coefficient 

e amplitude of the circular limit cycle 

F coulomb spline force vector 

I p equivalent polar moment of inertia 

k y , z equivalent stiffness in y,z direction 

l ll2 length of shaft 1,2 

m equivalent mass of the system 

t time, nondimensionalized using w y 

T torque transmitted by the coupling 

v,w response in the y,z direction (nondimensionalized relative to ) 

a gyroscopic parameter = I p /m if 

a 0 attitude angle of the elliptical response 
1 orthotropy parameter = w z / w y 

< viscous damping ratio 

x nondimensional nonlinearity parameter (16 splines) 

= 0.64/iT (l+ij/igj/mwya Jj* 

it coefficient of friction 

Q nondimensional spin speed of the rotor = Q/w y 
w nondimensional whirl speed of the rotor = w/w y 

w y A y /m 

w z ATT® 
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INTRODUCTION 


A spline coupling connecting two rotors is often used in the 
turbine engine industry. When it is not sufficiently lubricated, the 
resulting coulomb friction can produce a driving force that sustains an 
asynchronous whirling motion. 

A simple rigid rotor system coupled by such a spline was analyzed 
by Williams and Trent (ref. 1) using simulations on an analog computer. 
Morton (ref. 2) considered the effect of coulomb friction in drive 
couplings on the response of a flexible shaft and concluded that 
misalignments had the effect of stabilizing such a system. Marmol, 
Smalley and Tecza (ref. 3) modelled a coulomb spline by equivalent 
viscous friction and comparisons of their analytical results to those 
from a test rig were satisfactory. 

The model to be analyzed in this paper is very similar to the one 
in reference 1, and contains only two degrees of freedom. The supports 
have, in general, orthotropic stiffness and the model includes the 
gyroscopic effect. 

The emphasis here is on obtaining an improved understanding of the 
underlying phenomena rather than a detailed quantitative analysis of a 
large order system. It is for this reason that analytical and 
numerical-analytical methods have been applied to this problem (in 
addition to numerical simulation) so that a better insight could be 
gained. 

MATHEMATICAL MODEL 


A system consisting of two rigid rotors coupled by a coulomb 
spline is considered (fig. 1). They are supported on a mount with 
linear isotropic viscous damping and (in general) anisotropic 
stiffness. This is essentially a two degree-of-f reedom model and the 
equations of motion are found to be: 



The above equations are in nondimensional form and the deriva- 
tives are with respect to a nondimensional time. The spin speed of the 
rotor is denoted by Q, a denotes a gyroscopic parameter, $ the damping 
ratio and i represents an orthotropy parameter (f - 1: isotropic 

support). An equivalent force vector, F arises due to the friction 
excitation in the coulomb spline and is given by: 


F 


X 

^ — 1/2 
[ (w - Qv) 2 + ( v + Qw) 2 3 



( 2 ) 


where, x is a nondimensional parameter denoting the strength of the 
nonlinearity. 
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The system is thus modelled by a set of coupled, second order, 
nonlinear, autonomous equations ( eqn . 1) and a closed form analytical 

solution is difficult, if not impossible. 


SOLUTION: ISOTROPIC CASE 


When the support is isotropic (jT= 1), an exact solution can be 
obtained by assuming a circular response of amplitude e and a 
frequency w. 

i . e. , v = e cos wt 

w = e s in XJt 


(3) 


Substitution of this assumed solution form into equations (1) and 
(2) yields the following results: 


Q < w 


No Limit Cycle Possible 


Q > w 


Limit Cycle Ex is t s 


with w 
e 


Qa/2 + /l + (Qo/2 ) 2 
X/2«w 


(4) 

(5) 


It is interesting to note that the frequency of the limit cycle is 
the same as the undamped natural frequency of the linearized system and 
is independent of the strength of the nonlinearity (X) and the linear 
viscous damping ratio ($). The amplitude of the resulting circular 
limit cycle is directly proportional to X and reduces with increasing 
spin speed Q. Figure 2 shows the frequency spectrum of the numerically 
simulated response at a nondimensional spin speed of 3. Clearly, a 
single peak occurs at a frequency of 1,56 that agrees with the 
analytical prediction (eqn. 4). 

The onset of instability speed, at which a bifurcation into a 
limit cycle takes place, can be found to be 

Qo i = • . ( 6 ) 


Similar results were obtained in reference 1, but without the 
inclusion of the gyroscopic effect. 

/ — *• 

Figure 3 shows the force balance diagram for the isotropic system 
at a spin speed above the onset of instability. From energy considera- 
tions, it can also be shown that when the spin speed Q is less then the 
whirl speed w, the coulomb spline dissipates energy and hence cannot 
sustain a limit cycle. When Q > w however, the spline gives rise to 
equivalent negative damping. Then, the amplitude of the limit cycle 
can be determined by equating the energy dissipated by the linear 
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viscous damping to the energy added (to the system) by the coulomb 
spline forces. 

The stability of the limit cycle can be investigated by perturbing 
the limit cycle and writing linear variational equations for the 
perturbed motion. A detailed stability analysis indicates that the 
circular limit cycle is stable for all possible values of the 
parameters . 


SOLUTION: ORTHOTROPIC CASE 

When the support is not isotropic (Y * 1), the solutions are no 
longer amenable to an exact analysis. Hence, two approximate methods 
have been used and they are outlined below. 

Trigonometric Collocation Method (TCM) 

This is a numerical-analytic method that has been used quite 
widely for the solution of nonlinear boundary value problems. Ronto 
(ref. 4) has recently developed and formalized the method for applica- 
tion to periodic solutions and the basic philosophy is outlined below, 
omitting all the details and proofs. Given a set of nonlinear ODE’s 

x = 7 (x) , (7) 

a periodic solution x 0 (t) is sought in the form: 

. u ^ 

x 0 (t) = a 0 + ^ ( a„ cos nwt + b n sin nwt) (8) 

n = 1 

where the ith typical coefficients ¥ t and b t and the frequency w are 
unknown. This assumed solution is substituted into the differential 
equations and the resultant errors are forced to be zero at a 
prescribed number of time points (at least 2N+1). This leads to a 
sufficient number of nonlinear algebraic equations in the unknown 
coefficients and the frequency and can be solved using various 
numerical methods. 

Though this method was developed mainly for the solution of 
nonautonomous equations in reference 4, it has been successfully 
developed and applied for autonomous systems as well. 

Application of the TCM to the orthotropic system reveals some 
interesting results. The response is nearly elliptical (as can be 
expected) but als o has several higher freq uency compon ents. These 
higher harmonics have very low amplitudes and hence are not easily 
observed by numerical simulation. 

The effect of the higher harmonics is best seen in a semi-log plot 
of the frequency spectrum (fig. 4) that has_been generated by numerical 
simulation. The nondiraensional spin speed Q is equal to 3 and peaks 
are seen to occur at w, 3w, 5w, etc., where the fundamental frequency w 
is „ 2.92. This may be contrasted with the frequency spectrum for the 



isotropic case (fig. 2). The fundamental frequency of the limit cycle 
is seen to be independent of the strength of the nonlinearity (X) and 
the damping (5). The higher harmonic amplitudes are seen to increase 
with an increase in orthotropy. 

Numerical simulation has been carried out for a range of para- 
meters and the results compare very well with those obtained by the 
TCM. 

Analytical Approximation (AM) 

In order to obtain some idea about how the response variables vary 
with the system parameters, analytical expressions that approximate the 
response are quite useful. To derive these, since the higher frequency 
components have low amplitudes, the solution is assumed in the form of 
a single frequency elliptical component. As the system is autonomous, 
the time origin of the response is arbitrary. Hence, 

v = a a cos ( w t + a 0 ) + a d cos(wt - « 0 ) ^ 

w = a m sin ( wt + « 0 ) - a d sin(wt - « 0 ) 

Substituting this assumed form of the solution into the equations 
of motion, equation 1, and after substantial algebraic manipulations 
and retention of only first order approximations, an approximate 
expression for the fundamental frequency of the limit cycle is seen to 
be : 




w2 « A {1 + /I - (//A)* } 

(10) 

where : 


A = (1 + r* + «2Q2)/2 

(ID 


and, 

b/a * (w 2 — l)/aQw 

(12) 



b « \/2«w 

(13) 


The approximate expression for the fundamental frequency is the 
same as the undamped natural frequency of the linearized orthotropic 
system. The general variation of the frequency and amplitude is very 
similar to the isotropic case. However, in contrast to the isotropic 
case, the limit cycle appears to exist at all spin speeds. This is 
consistent with the simulation results of Williams and Trent (ref. 1). 

Table 1 shows a comparison between the analytical predictions and 
the values obtained by the TCM for several values o f f (o rthotropy) . 
The frequency prediction is correct to within 10% for f up to 3.0 and 
the amplitude prediction has less than 10% error for i up to 1.6. 

The stability of the limit cycle in the orthotropic case has not 
yet been analyzed and further investigations will be reported 
elsewhere . 
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CONCLUSION 


A two degree-of-freedom model of rigid motors coupled by a coulomb 
spline is analyzed. Two different cases are considered. 

When the support is isotropic, a bifurcation into a stable, 
circular limit cycle takes place at a spin speed equal to the forward 
undamped critical speed. The frequency of the limit cycle is the same 
as the undamped natural frequency of the system and is independent of 
the strength of the nonlinearity. 

When the support is orthotropic, the limit cycle has many 
frequency components but is still nearly elliptical. The undamped 
natural frequency of the system, however, approximates the fundamental 
whirl frequency of the response to a very good accuracy. The higher 
harmonic amplitudes increase with an increase in orthotropy. The 
stability of the response in the orthotropic case is of considerable 
interest and remains to be investigated. 

Similar qualitative phenomena can be expected to occur in large 
flexible rotor-bearing systems with a spline coupling. Work is 
continuing along these lines. 
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TABLE 1 . - COMPARISON BETWEEN ANALYTICAL AND TCM RESULTS 



( Q = 3. 

0, x = 0 

.01867, 

« = 0.1, 

a = 0.3) 


1 

“am 

“tcm 

(b/a) A „ 

(b/a) TCM 

^AM X ^® 2 

^TCM X ^ 2 

1.0 

1.5466 

1.5466 

1.0000 

1.0000 

6.0348 

6.0348 

1.1 

1.5470 

1.570 

1.0007 

1.0007 

5.8511 

6.0169 

1.2 

1.6495 

1.6490 

1.1591 

1.1513 

5.6584 

5.8875 

1.3 

1.7093 

1.7081 

1.2491 

1.2348 

5.4605 

5.6251 

1.5 

1.8433 

1.8386 

1.4453 

1.3929 

5.0633 

4.6291 

3.0 

3.1466 

2.8387 

3.1431 

1.1148 

2.9662 

0.4871 


y 
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